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Abstract 

We consider the influence on the surface critical behaviour of a quantum Ising chain 
of quenched random surface perturbations decaying as a power of the distance 
from the surface (random Hilhorst-van Leeuwen models). We study, analytically 
and numerically, the multiscaling behaviour of the surface magnetization and the 
surface energy density in the case of marginal perturbations. 
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1 Introduction 



In semi-infinite two-dimensional (2D) Ising models, extended perturbations of 
the coupling constants of the form AKi = Al~ K , decaying as a power of the 
distance I from a free surface, were first studied by Hilhorst and van Leeuwen 
(HvL) at the begining of the 1980s [1]. Since then they have been the subject of 
continuous interest [2-15]. Similar extended perturbations associated with line 
defects in the bulk [16-18] or point defects [19-21] have also been considered 
(see Ref. [22] for a review). 

Such extended perturbations are known to modify the local critical behaviour, 
for an arbitrarily small value of the perturbation amplitude A, as soon as the 
decay exponent k is sufficiently small. Analytical results for a surface defect in 
the 2D Ising model [1,6] show that with a temperature-like perturbation (the 
case considered later) marginal behaviour is obtained when k — 1. Then the 
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surface exponents vary continuously with the perturbation amplitude A and 
the surface transition becomes first order when A is greater than a critical 
value A c . When k < 1, the surface properties display essential singularities 
when A < and the transition is first order when A > 0. 

This rich critical behaviour has been explained independently by Cordery and 
Burkhardt [2,3] through scaling considerations. Under a change of the length 
scale /' = l/b, the perturbation transforms as 



where v is the bulk correlation length exponent. Comparing both sides of the 
last equation, one obtains the following scaling behaviour for the perturbation 
amplitude: 



Thus, for the 2D Ising model, k c — 1/v — 1 corresponds to the critical value of 
the decay exponent below which the extended perturbation becomes relevant. 

An interesting aspect of these smoothly inhomogeneous models is that in 
2D their correlation functions are covariant under conformal transformations 
when the perturbation is marginal. Thus gap-exponent relations are satis- 
fied [12,14,18-20] and conformal profiles are obtained in the strip geome- 
try [15], provided the form of the perturbation is also properly transformed as 
first noticed in Ref. [12]. 

Recently, aperiodic [23] and random [24,25] versions of the HvL model have 
been studied for the quantum Ising chain in a transverse field, which corre- 
sponds to the extreme anisotropic limit [26] of the 2D classical Ising model 
and belongs to the same universality class. 

In Ref. [24] random extended perturbations decaying towards a pure bulk 
system were considered with two types of quenched randomness. 

In model #1, the random couplings take the following form: 



[A*i] 



' = — = h 1/u — 



(1) 



A' = b- K+1/u A . 



(2) 





with probability 



1/2 
1/2 



(3) 



Thus the average perturbation (J; — J) vanishes. 
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In model #2 

J t = J[l + /,a 2 ] , 



with probability 



1 - Z~ K 



(4) 



Then the average coupling (Ji) = J(l + a 2 ^ K ) has the same form as in the 
HvL model. 

In both models the transverse fields are kept constant, hi = h. 

A perturbation expansion of the average free energy and scaling considerations 
allowed us to obtain a relevance criterion for both types of perturbations [24]. 
With model $4, the first-order correction to the average free energy vanishes 
and the second-order correction leeds to the following scaling behaviour for 
the amplitude a± of a thermal perturbation: 



With model #2, the scaling behaviour of the perturbation amplitude a 2 is 
governed by the non-vanishing first-order correction to the average free energy 
so that 

a ' 2 = b^- K a 2 , (6) 



as for the non-random HvL model in (2). According to (5) and (6) marginal 
behaviour is obtained for k — 1/2 with model #1 and k—1 with model #2. 

Analytical expressions for the finite-size scaling behaviour at criticality of the 
average and typical surface magnetization were obtained in Ref. [24], in the 
case of marginal and relevant perturbations, in agreement with the results 
of numerical simulations. We gave also some conjectured expressions for the 
surface energy exponents in the marginal case. 

In the present work, we extend our previous study by considering the finite-size 
scaling behaviour of the average of the moments of the surface magnetization 
(m|(L)) and of the singular part of the surface energy (ej(L)) in the marginal 
situation. In Section 2, we briefly recall the techniques used to diagonalize the 
corresponding quantum Hamiltonians and to obtain the scaling behaviour of 
the surface magnetization and surface energy. Analytical expressions for the 
average of the moments of the surface magnetization and the surface energy are 
deduced from scaling considerations in Sections 3. The results are compared 
to numerical finite-size scaling data in Section 4 and their range of validity 
discussed in Section 5. 



3 



2 Free fermions techniques 



2. 1 Diagonal Hamiltonian 



We consider the semi-infinite transverse-field Ising model (TIM) with Hamil- 
tonian 



n = -lJ2(J l afaU 1 + h l af), (7) 
i=i 



where of and of are Pauli matrices, J\ is the random coupling and h\ the 
transverse field. 



The TIM Hamiltonian can be diagonalized using standard methods [27,28]. 
A Jordan-Wigner transformation [29] changes the original Hamiltonian (7) 
into a quadratic form in fermions which is diagonalized through a canonical 
transformation leading to 



n = Y,e q {r,\r, q - \) (8) 

0=1 



for a chain with length L. The ry (r] q ) are fermion creation (annihilation) op- 
erators and the non-negative excitation energies e q satisfy the set of equations 



e q ij) q (l) = -h l <l> q (l)-J l <i> q (l + i), 

e q <f>S) = -Ji-iM ~ !) - Witt) (9) 



with the boundary conditions Jo = Jl = 0. Introducing 2L-dimensional vec- 
tors V q with components 



V q (2l-l) = -<f> q (l), V q (2l)=$S) ( 10 ) 
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the relations (9) lead to an eigenvalue problem for the following tridiagonal 
matrix: 

' fn 
hi Ji 
t _ Ji h 2 

Ji-i h L 
v h L 

When T is squared, odd and even components of V q decouple and one obtains 
two separate eigenvalue problems for (ft and The relations (9) are invariant 
under the transformation <p q — > — <p q , e q — > —e q so that, changing <p q into 
— <fi q in V q , one obtains the eigenvector of T with eigenvalue — e q . Thus the 
physically relevant information is contained in that part of the spectrum of T 
with e q > 0. 



2.2 Surface magnetization and surface energy 

The surface critical properties can be obtained through finite-size scaling at 
criticality, i.e., by working on a finite system with length L and J = h = 1. 

The imaginary time spin-spin autocorrelation function is given by 

QT(r) = ((K(tK(0)|0> = £ \(n\af\0)\"eM-r(E n - E Q )} , (12) 



where |0) and \n) are the ground state and the nth excited state of Ji, E and 
E n the corresponding eigenvalues. With fixed boundary conditions, o z L = ±1, 
at I — L, the ground state is degenerate and asymptotically, one obtains 
lim^oo Qf a {T) = mf where 

mi = (a\af\0) (13) 

is the local magnetization. It is given by an off-diagonal matrix element in- 
volving the first excited state with one fermionic excitation, \a) = f]\\0), which 
is degenerate with the ground state. For the surface spin at / = 1 a simple 



(11) 
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expression is obtained [5]: 



m s (L) =m l = 



1=1 k=l \ J kJ 



-1/2 



(14) 



This expression can be rewritten under the following form [24]: 

L-l 



m s (L)=m^(L)n 



Ji 



Li W ' 



(15) 



where 



ff^(L) 



L-l L-l 

e n 

3=1 k=L-j 



-1/2 



(16) 



is the surface magnetization at / = L — 1 on the dual chain with transverse 
fields Ji, couplings hi (I = 0, L — 1) and fixed boundary conditions Jo = at 
1 = 0. 

The scaling dimension of the surface energy density can be obtained by consid- 
ering the finite-size behaviour with free boundary conditions of the off-diagonal 
matrix element [30] 



e B (L) = (e|(7f|O> = (e 2 -e 1 )^(l)0 2 (l), 



(17) 



where \e) = TjlrjllO) is the lowest eigenstate leading to a non- vanishing ma- 
trix element. This matrix element enters into the expression of the connected 
energy-energy surface autocorrelation function Gf £ (r) and scales like the sin- 
gular part of the surface energy density. 

For the critical HvL model with a marginal decay of the perturbation, hi — 1, 
Ji = 1 + a/l, the product in (15) behaves asymptotically as L a . When L ^> 1, 
the couplings on the right side of the chain are asymptotically unperturbed, 
thus one expects that rfi^(L) ~ L -1 / 2 as for the homogeneous Ising model at 
the ordinary transition. This leads to the scaling behaviour: 

m s (L)~L"<, x s m = l-a, a<\. (18) 



This expression cannot be valid beyond a = 1/2 where the surface remains 
ordered at the bulk critical point [6,5]. The surface magnetization then displays 
a first order transition and = when a > a c = 1/2. Thus our assumption 
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concerning the behaviour of fn^(L) is wrong when a > 1/2 and then one must 
have TrT^(L) ~ L~ a . 

For the surface energy in Eq. (17), when a < 1/2 both 0i(l) and 02(1) scale 
as m s (L) ~ L -3 ^. The excitations ei and e 2 scale as L^ 1 leading to 

e s (L) ~ L-'^(i) ~ a* =2(1 -a), a<§, (19) 

in the second-order regime. When a > 1/2, the onset of surface order leads 
to a more complicated scaling behaviour. The first gap e\ is associated with a 
localized state which reflects the surface behaviour. It scales as L~ 2a while e 2 ~ 
L^ 1 as usual. We have 0i(l) = rn s (L) ~ L° but 02(1) behaves anomalously as 
L~ a+l / 2 . Finally, x% — a + 1/2 in the first order regime, a > 1/2. 



3 Calculation of the moments 

In this section we evaluate the moments of the surface magnetization and the 
surface energy density at criticality on a chain with length L in the regime of 
second-order surface transition. We restrict our study to the case of marginal 
perturbations, i.e., with k — 1/2 for model #1 and k — 1 for model #2. 



3.1 Surface magnetization 

We assume that the surface magnetization at I = L — 1 on the dual chain, 
rn^(L) in (15), is still scaling as L -1 / 2 for the random models so that 

k(l)> = n j?) = l-^ nVf ) , (20) 

1=1 1=1 

where we used the statistical independence of the couplings on different bonds. 
Using the distribution of the couplings given in (3) for model #1, with J = 1 
we obtain 

w^+i&y+K'-i&y^+^f (2i> 

The last expression gives the asymptotic behaviour when / ^> 1 which governs 
the critical behaviour. Taking the logarithm of the product in (20) and a 
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continuum approximation leads to 

m M) - /m (i + ) di , Mi^M^i. (22) 

Thus the gth moment of the surface magnetization scale as 

(m q s {L)) ~ ~ £-(<//2)[i-(9-i)a?] (23) 

and the surface magnetization displays multiscaling with 

<(?) = |[l-(9-l)«?] (model #1). (24) 

Similarly for model #2, using (4) with J = 1 and k = 1, we obtain 

W .) = i(l + n2 )< + 1 -l = 1+ ( 1 + a f- 1 . (25) 



Since 

L ^ ' r (l + a 2 )*-l 



i=i i 



i+ 



i 



d/~ [(l + a 2 )"-l]lnL, (26) 



the qth moment of the surface magnetization scales as 

K(L)) ~ L-(*/2)+(i+a2)«-i , (27) 



so that 



4.(9) = \~ {l + a2)q 1 (model #2) . (28) 



The exponents in Eqs. (24) and (28) are in agreement with the average and 
typical exponents obtained previously in Ref. [24] when q — 1 and q — > 0, 
respectively: 

KJ av = |, K] tyP = |(l + «i) (model #1) 

= I - 02 , K] tyP = I - Ml + 02) (model #2) 
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3.2 Surface energy density 



For the scaling behaviour of the moments of the surface energy density, one can 
start from the expression giving e s (L) in Eq. (19). The factor L~ x is associated 
with the scaling of the excitations and reflects the isotropy of the bulk critical 
behaviour. It should not be modified with a random surface perturbation, 
hence we have 

(ei(L)) ~ L-i{ml«{L)) . (30) 

The gth moment of the surface energy density depends on the moment of 
order 2q of the surface magnetization. 

For model #1, making use of Eq. (23), we obtain 

(e q s (L)) ~ L- qx * {q) ~ L- q[2 -( 2q - l)a * ] , (31) 

so that 

x s c (q) = 2-(2q- l)a\ (model #1) . (32) 

In the same way Eq. (27) leads to 

(e q (L)) ~ ^+(1+^-1 (33) 



for model #2 so that 



xl(q) = 2- (1 + fl2)2 ' (model #2) . (34) 



The average and typical exponents are then given by: 

(model SI) 

(35) 



[^] av = 2 - a\ , [^] typ = 2 + a\ (model #1) 



= 2 - a 2 (2 + a 2 ) , [x% yp = 2[1 - ln(l + a 2 )\ (model #2) 



The typical exponents in Eqs. (29) and (35) are in agreement with the values 
conjectured in Ref. [24]. They keep the same form as for the pure HvL model 
in Eqs. (18) and (19) with a replaced by an effective amplitude, —a\/2 for 
model SI and ln(l + 02) for model S2. These typical effective amplitudes are 
deduced from the relation In Ji(a e g) ~ (In J/) for I ^> 1. A similar conjecture 
proposed in Ref. [24] for the average exponents is actually verified only up to 
the first order in the perturbation amplitude. 
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4 Numerical results 



4-1 Surface magnetization 

The surface magnetization for a given configuration of the couplings {Jk} was 
obtained using Eq. (14) at criticality, i.e., with — 1 and J = 1 in (3) 
and (4). 

1.5 
1.0 

V 

0.5 



0.0 

0.0 0.2 0.4 0.6 0.8 1.0 

a 1 

Fig. 1. Exponent of the qth. moment of the surface magnetization as a function of the 
perturbation amplitude for model #1. The extrapolated exponents were deduced 
either from exact enumerations for chains with sizes L = 4-24 (open symbols) or 
from Monte Carlo samplings for chains with sizes L = 2 5 -2 14 (filled symbols). The 
lines correspond to the analytical expression in (24). 

The moments were averaged either by performing exact enumerations of the 
random configurations in the case of chains with lengths L = 4-24 or through 
Monte Carlo samplings (over 10 6 samples) for longer chains with lengths 
L = 2 5 -2 14 . The exponents x s m (q) were deduced from an extrapolation of two- 
point approximants using the BST algorithm [31]. The marginal exponents 
are shown in Fig. 1 for model #1 and in Fig. 2 for model #2. 

4-2 Surface energy density 



The surface energy density was obtained using the expression given in (17) 
where e±, e 2 , 0i(l) and 02(1) are obtained through the diagonalisation of a 
L x L matrix constructed from the elements of T 2 with odd indices. 
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-0.9 -0.6 -0.3 0.0 0.3 0.6 

a 2 

Fig. 2. As in Fig. 1 for model #2. The lines correspond to the analytical expression 
in (28). 




a 1 



Fig. 3. Exponent of the ^th moment of the surface energy density as a function of the 
perturbation amplitude for model #1. The extrapolated exponents were deduced 
either from exact enumerations for chains with sizes L = 4-24 (open symbols) or 
from Monte Carlo samplings for chains with sizes L = 2 2 -2 9 (filled symbols). The 
lines correspond to the analytical expression in (32). 

The average moments were deduced from exact enumerations for chains with 
lengths L = 4-24 and from Monte Carlo samplings (over 10 4 samples) for 
chains with lengths L = 2 2 -2 9 . The exponents x s c (q) were deduced from an 
extrapolation of two-point approximants using the BST algorithm in the case 
of exact enumerations. The approximants deduced from the Monte Carlo data 
were too noisy to use the BST extrapolation process. In this case the exponents 
were deduced from a non-linear fit of the following expression: 

ln(e«(L)> ~ In A - x s e (q) In L + BL^ . (36) 
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1.0 

-0.9 -0.6 -0.3 0.0 0.3 0.6 

a 2 

Fig. 4. As in Fig. 3 for model #2. The lines correspond to the analytical expression 
in (34). 

The continuously varying exponents are compared to the analytical expres- 
sions in Fig. 3 for model #1 and in Fig. 4 for model #2. 



5 Discussion 

1.0 
0.9 
0.8 

5 °- 7 

s 

'* 0.6 
0.5 
0.4 

-0.5 -0.2 0.0 0.2 0.5 0.8 1.0 

a. 

Fig. 5. Exponent of the qth moment of the dual surface magnetization at L as a 
function of the perturbation amplitude for model #1 (full symbols) and model #2 
(open symbols). The extrapolated exponents were deduced from exact enumerations 
for chains with sizes L = 4-24. The line corresponds to the behaviour at the ordinary 
surface transition where = 1/2. 

A good agreement between numerical and conjectured analytical results have 
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been obtained for the surface magnetization and the surface energy density, 
except in the vicinity of first-order surface transitions. 

Already for the pure HvL model a slow convergence to the exactly known 
behaviour was observed close to the first-order regime [13]. One may verify 
on Figs. 2 and 4 for model #2 that the Monte Carlo data, obtained on larger 
systems, are closer to the analytical results when 02 > 0. Nevertheless, like for 
the the pure system, some of the scaling assumptions leading to the analytical 
expressions for the surface exponents cannot be valid beyond some critical 
value of the perturbation amplitude since x s m (q) and x s e (q) have to remain 
non- negative. 

The dual magnetization in (15) was assumed to scale as L~ 1//2 , like for the 
unperturbed system at the ordinary transition. We performed a numerical 
study of the moments of Tn^(L) for the two models using Eq. (16), exact 
enumerations for chains with lengths L = 4-24 and BST extrapolations. The 
results are shown in Fig. 5. 

For model #1, the ordinary surface behaviour is obtained for any value of 
the perturbation amplitude which gives some confidence in the conjectured 
analytical expressions in the regime of second-order transition. For model #2, 
the ordinary surface behaviour is only obtained when a<i < .2. Beyond this 
value deviations from the ordinary surface behaviour are observed and the 
analytical expressions for the exponents are probably no longer valid. 
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